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ABSTRACT 

We consider the origin of the so-called S stars orbiting the supermassive black 
hole at the very center of the Galaxy. These are usually assumed to be massive 
main-sequence stars. We argue instead that they are the remnants of low-to- 
intermediate mass red giants which have been scattered on to near-radial orbits 
and tidally stripped as they approach the central black hole. Such stars retain 
only low-mass envelopes and thus have high effective temperatures. Our picture 
simultaneously explains why S stars have tightly-bound orbits, and the observed 
depletion of red giants in the very center of the Galaxy. 

Subject headings: Galaxy: center — Galaxy: kinematics and dynamics - Galaxy: 
nucleus 


1. Introduction 

Proper motion observations (Eckert & Genzel 1997, Ghez et al 1998) have revealed in 
exquisite detail the motions of stars around the black hole in the center of the Galaxy. The 
closest known stars have extremely eccentric orbits (for example S2 has an eccentricity of 
0.87), pericenters of only a few hundred AU, and orbital periods in the range 15 - 300 yr 
(Schodel et al 2003, Eisenhauer et al 2005). These “S stars” are all rather blue (T e g ~ 
30000 K. These properties lead to the “birth problem” (Ghez et al 2003). The S stars 
orbit in a region where tidal forces from the central supermassive black hole prevent star 
formation in situ, although several authors have examined whether they were produced after 
a collision of dense molecular clouds or in an AGN-likc disc (e.g. Genzel et al. 2003; 
Levin & Beloborodov 2003; Milosavljevic & Loeb 2004). A plausible alternative picture is 
that the S stars result from the sinking of massive stellar clusters towards the black hole 
by dynamical friction. However tides disrupt such clusters at distances ~ 1 pc from the 
black hole, much further out than the region of the observed S stars (Gerhard 2001). To 
supply the S stars requires scattering into near-radial orbits by gravitational interactions 
with other stars. However this process occurs on the local relaxation timescale, which would 
considerably exceed the main-sequence lifetime (around 10 Myr) of stars of the observed 
temperatures. 
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One way out of this problem is to increase the stability of the clusters against tides. 
There have been several suggestions that stellar clusters could contain intermediate-mass 
black holes, making them less vulnerable to tidal shredding (Hansen & Milosavljevic 2003, 
McMillan & Portegies Zwart 2003, Kim et al 2004). We do not adopt this idea in this Letter. 

These closest stars have another intriguing property: there appear to be no red giants 
among them, despite the fact that these would be brighter than main-sequence stars in 
the infrared K band where the observations are made. We call this the red-giant depletion 
problem. It has been suggested that such stars are destroyed in stellar collisions (e.g. Genzel 
et al. 1996; Alexander 1999) although collisions have been shown to be inefficient at removing 
most of the red-giant envelope (Bailey & Davies 1999). 

A third problem arises from a selection effect for identifying stars moving very close to 
the black hole pointed out by Schodcl et al (2003). Namely, stars with pericenters similar to 
the S stars but much longer orbital periods would spend almost all of their lifetimes far from 
pericenter moving very slowly, and thus not turn up in proper motion observations. In other 
words, only tightly bound orbits are easily identifiable, and a valid picture of the origin of 
the S stars must explain these orbits too. We call this the orbit problem. 

In this Letter we suggest a simultaneous solution to the birth, red-giant depletion and 
orbit problems. 


2. Violent Mass Loss 

Table 1 gives the known pericenter distances p and eccentricities for S stars. The 
instantaneous tidal lobe of a star at pericenter has size 


R t ^ 0.5 
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Mi 
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where M*, Mbh are the masses of the star and central black hole. Table 1 gives the values of 
Rt corresponding to the observed S star orbits, taking values of 3M 0 and 3.6 x 1O 6 M 0 for 
M*,Mbh- hi four cases (including those with the most accurate orbital determinations) we 
note that these are of order a few 1OOR 0 . This is comparable in size to a low-to-moderate 
mass star on the AGB, with core masses ~ 0.6 — 0.8 M 0 . This can be seen in Fig. 1, where 
we plot the evolutionary tracks for solar metallicity stars of masses 1, 2, 3, and 4 M 0 using 
the SSE code of Hurley, Pols & Tout (2000). 


We can now ask what would happen to an AGB star scattered on to an orbit passing 
close enough to the black hole that its radius exceeds its tidal radius. Tidal disruption of 
stars in the center of our own and other galaxies has been discussed in some detail (e.g. 
Alexander & Livio 2001; Di Stefano et al 2001; Ghez et al 2003). At its first pericenter 
passage, the scattered star loses most of its envelope in a violent fashion, leaving only the 
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carbon-oxygen core, the shell-burning layers, and a tenuous envelope. Observations of the 
secondary star of the X-ray binary Cyg X-2 (Casares et ah, 1998; Orosz & Kuulkers 1999) 
show that it still has hydrogen in its envelope, and that it has a luminosity L ~ 150 L 0 
despite a mass of only ~ O.8M 0 . These properties agree with binary evolution models where 
the secondary descends from a 3.5 M 0 star that has lost most of its mass quite rapidly (King 
& Ritter, 1999, Kolb et al, 2000, Podsiadlowski et al 2002). These models suggest that the 
envelope mass is less than 0.1 M 0 on top of a 0.7 M 0 core, but the star closely resembles a 
3.5 M 0 main-sequence star with L ~ 150 L 0 . We expect that an AGB star that encounters 
the black hole will be left with a similar envelope mass. The star retains the same nuclear 
luminosity but is much smaller and bluer. Stripped AGB stars have luminosities around 10 4 
L 0 and temperatures > 10 4 K, similar to B0 main-sequence stars. This is consistent with 
the idea that the S stars were all red giants which were violently stripped of their envelopes 
near pericenter. The nuclear lifetimes of these stripped stars would be a few Myr, since they 
have < O.1M 0 of envelope available for nuclear burning (see the Discussion). Stars must be 
captured and stripped at a rate ~ 10~ 6 yr _1 to prevent the current population of S stars 
decaying on this timescale. Since the tidal interaction is effectively pointlike, the percenters 
of the stripped stars remain unchanged by it. As seen from the Table, this kind of evolution 
is possible for the current orbits of at least four of the S stars. We discuss the other two 
stars later (Section 3). 

Stripping stars solves the birth problem: if the S stars are not main-sequence stars 
with masses ~ 10 — 15M 0 , but instead stripped giants of orginal (main-sequence) masses 
~ 1 — 4M 0 , their ages are greater than 300 - 3000 Myr, allowing ample time for them to be 
scattered from star clusters further out into orbits with the current small pericenter distances. 
The violent tidal mass loss we have invoked leaves the pericenter distances unchanged, but 
produces the observed tightly bound orbits, as orbital energy is used to remove the stellar 
envelope. 

To see this we compute the energy loss required to take a star from a near-parabolic 
orbit to a tightly bound one like those seen, and compare it to the energy needed to strip a 
red giant down to a core of mass M c and radius R c . The loss of orbital binding energy is 
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where a is the semi-major axis after the mass loss. To strip the star down to the core plus 
a burning shell and a very tenuous envelope requires an energy 
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where M e is the initial envelope mass. Thus stripping requires 
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The resulting values of a compare well with those measured for the S stars listed in Table 
1, and imply periods in the observed range < 100 yr. The observed percenters p = a(l — e) 


require eccentricities 


e ~ 1 — 2 


pM e 
R c M b h 


~ 0.9 


( 5 ) 


The idea of violent tidal mass loss thus explains the current tight orbits of most of the 
S stars as well as solving the birth problem. We can extend this to a simultaneous solution 
of the red-giant depletion problem by noting our conclusion at the end of the Introduction 
above that only tightly-bound orbits are easily identifiable. Red giants on orbits such that 
they do not fill their tidal lobes near pericenter have periods of many centuries and are very 
unlikely to be identified as high proper motion objects. Hence the only stars potentially 
identifiable as Galactic Center red giants lose their envelopes and turn into S stars instead. 


3. Scattering of Stars 


We noted above that two of the S stars (SI and S13) would not have filled their tidal 
lobes as red giants at the current pericenters, despite having similar semi-major axes and 
thus binding energies. These two orbits have significantly smaller eccentricities, implying 
higher angular momenta. This strongly suggests that they have been scattered from their 
original post-interaction orbits by the addition of angular momentum near apocenter. For 
this to work, the mean gravitational relaxation time should not be much greater than the 
nuclear lifetimes ~ few x 10 6 yr of these stripped stars. As a guide we use the estimate 

+ n o/i c’’ 3 

W = °- U G 2 Mp\nA 


(Binney & Tremaine, 1987, eq. (8-71), where a is the velocity dispersion, M the mass of a 
star, and p the stellar mass density. Strictly, this formula needs some modification for the 
presence of the central black hole, but this should not change things by large factors.) If stars 
have been fed into the central regions for a significant fraction of the Galaxy’s lifetime there 
could currently be ~ 10 -6 x 10 10 = 10 4 remnant stars (mostly white dwarfs, but with a few 
neutron stars and black holes) orbiting within the central 0.05 pc, giving p ~ 10 8 M e pc -3 . 
The remnants spend almost all their lifetimes near apoastron, so we take cr of order the 
orbital velocities ~ 150 km s _1 there. This gives f re i ax ~ 10' yr. This supports the idea 
that the orbits of SI and S13 have been scattered once from their original post-interaction 
forms. We note that these orbits are not relaxed, and still fairly close to the original ones, 
as expected if the number of scatterings is low. 

Stellar-mass black holes produced in massive supernovae may settle into the central 
parsec by dynamical friction (Morris 1993; Miralda-Escude & Gould 2000). In doing so they 
may scatter the white dwarfs out of the central parsec. The resulting central density will be 



similar to that calculated above (see also Alexander & Livio 2004), but the average mass of 
the stars would be a factor of ten higher, which would shorten the relaxation time. 


4. Discussion 

We have shown that the S stars can result from violent stripping of red giants near peri- 
center, making them into rather blue objects. Their subsequent evolution is straightforward. 
Their nuclear luminosities are essentially unaffected by the loss of much of the envelope, so 
they evolve at constant luminosity until the remaining envelope mass is burnt. This takes of 
order a few xlO 6 yr. Their radii shrink and their effective surface temperatures T increase 
during this phase. The constant total luminosity fixes AT 4 (A = area), while the K-band 
flux varies as ~ AT, so they become fainter in the K-band as ~ T~ 3 . This is probably 
the reason why all the observed S-stars appear to have very similar temperatures. A star 
stripped on the RG branch (i.e. pre-helium ignition) has a luminosity ~ 10 2 L 0 , and with 
an effective temperature T e g ~ 10 4 K would not be visible in IR surveys. We expect the 
S stars to remain observable for some fraction of the nuclear timescale rsj 10 6 yr. 

To maintain a steady state, S stars must be created at the same rate ~5x 10 -6 yr^ 1 
(assuming of there are of order 5 at any moment). This fixes the rate at which stars must be 
scattered into near-parabolic orbits with the very small pericenters observed for the S stars. 
The total stellar tidal disruption rate by the central black hole is ~ 5 x 10 -4 yr -1 , scaling 
as Rl /4 (Wang & Merritt 2004). Considering the evolutionary tracks plotted in Fig. 1, we 
find that about 1 % of stripped 1-4 M 0 stars are stripped while on the AGB. Hence the 
tidal disruption rate of AGB stars agrees well with the formation rate of the S stars in our 
picture. 

The tidal disruption rate derived above implies that gas is injected into orbit about the 
black hole at a rate ~ 2.5 x 1O _4 M 0 yr -1 . This is significantly higher than estimates of the 
black hole accretion rate (Baganoff et al 2003). However this is not suprising, as the gas still 
has to lose almost all its angular momentum in order to accrete. It is likely that accretion 
is a highly time-dependent phenomenon, with the Eddington limit playing a significant role 
in determining the growth rate of the black hole mass. 

We note finally that our picture implies the presence of a kind of “Oort cloud” of stellar 
remnants consisting mainly of white dwarfs around the central black hole. The dynamical 
evolution of this system may be very interesting. Moreover the stars in this system must 
interact with the orbiting gas discussed in the last paragraph. This may be the origin of 
the observed X-ray flares in a similar manner to that suggested by Nayakshin & Sunyaev 
(2003). We might expect the gas to virialise with the stellar population, implying X-ray 
temperatures (~ 10 8 K) near the black hole. 

In a quite general way, our picture suggests that the stellar populations near the Galac- 
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tic Center should show orbital, positional, and surface characteristics determined by their 
interaction with the central black hole. We will return to these questions in future papers. 
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Table 1. 


Tidal radii for measured orbits. 


Name 

p (mpc) 

e 

Rt (r©) 

p (mpc) 

e 

Rt (Rq) 

S2 

0.6 

0.87 

114 

0.6 

0.88 

111 

S12 

2.2 

0.73 

424 

1.1 

0.90 

206 

S14 

0.4 

0.97 

75 

0.5 

0.94 

101 

SI 

8.0 

0.62 

1541 

10.1 

0.36 

1941 

S8 

0.6 

0.98 

116 

0.9 

0.93 

176 

S13 

6.0 

0.47 

1156 

5.0 

0.40 

963 


Note. — The 2nd-4th columns are from the orbital fits of Schodel 
et al (2003), while the 5th-7th columns are from the fits of Eisen- 
hauer et al (2005). Note that the orbits for S14, SI, S8, and S13 
were all poorly constrained in Schbdel et al. 







